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Abstract. Inspired by Anderson et al. [J. R. Soc. Interface, 2014] we 
study the long-term behavior of discrete chemical reaction networks 
(CRNs). In particular, using techniques from both Petri net theory and 
CRN theory, we provide a powerful sufficient condition for a structurally- 
bounded CRN to have the property that none of the non-terminal reac¬ 
tions can fire for all its recurrent configurations. We compare this result 
and its proof with a related result of Anderson et al. and show its con¬ 
sequences for the case of CRNs with deficiency one. 


1 Introduction 

Chemical reaction network (CRN) theory studies the behavior of chemical sys¬ 
tems. Traditionally, the primary focus is on continuous CRNs, where mass action 
kinetics is assumed, see, e.g., mm- In this setting a state is determined by the 
concentration of each species and the system evolves through ordinary differen¬ 
tial equations. However, in scenarios where the number of molecules is small one 
needs to resort to discrete CRNs. In a discrete CRN a state (also called config¬ 
uration) is determined by the counts of each species, and one often associates a 
probability to each reaction. In this paper we consider only discrete CRNs, and 
so, from now on, by CRN we will always mean a discrete CRN. 

A CRN essentially consists of a finite set of reactions such as A + B —?> 
2B, which means that during this reaction one molecule of species A and one 
molecule of species B are consumed and as a result two molecules of species B 
are produced. We may depict a CRN as a graph, the reaction graph, where the 
vertices are the left-hand and right-hand sides of reactions and the edges are the 
reactions, see Figure [T] for an example. We focus in this paper on the long-term 
behavior of CRNs for which the number of molecules cannot grow unboundedly. 
For such CRNs, called structurally-bounded CRNs, each configuration eventually 
reaches a configuration c such that c is reachable from any configuration c' 
reachable from c (i.e., we can always go back to c). Such configurations are 
called recurrent. The CRN N of Figure [T] is structurally-bounded. 

Now, let us consider the CRN N' obtained from N by replacing every vertex 
by one molecule of a distinct species Xi, see Figure [5J We easily observe that for 
N', the recurrent conhgurations are exactly those without molecules of species 
Xi or A 5 . In other words, the reactions Pi and P 5 cannot fire for any recurrent 
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A + B 2B —^ B + C B —^ A -- D 



2C 


Fig. 1. The reaction graph of a CRN N. 



X4 


Fig. 2. The reaction graph of the CRN N' obtained from N by introducing a distinct 
species Xi for each vertex. 


configuration of N'. Notice that the reaction graph of N' has two strongly- 
connected components without outgoing edges: one having the vertices X2, X3, 
and X4 and one having the vertices Xq and Xj. The reactions outside these 
two strongly-connected components are called non-terminal. Thus N' has the 
property that none of the non-terminal reactions can fire for all its recurrent 
configurations. But what about the original CRN iV? The dynamics of N are 
clearly more involved since we can go, for example, from configuration A + B 
back to A + B hy firing reaction ai followed by firing reaction . 

The main result of this paper, cf. Theorem [l|, is a sufficient condition for a 
structurally-bounded CRN to have the property that none of the non-terminal 
reactions can fire for all its recurrent configurations (we recall the notion of non¬ 
terminal reaction in Section]^). Those CRNs have relatively simple long-term 
behavior. The sufficient condition of Theorem [3 (when formulated in terms of so- 
called T-invariants in Corollary @) is structural/syntactical and can be checked 
for many CRNs in a computationally-efficient way. Various non-trivial CRNs 
from the literature satisfy the sufficient condition of Theorem [l| (see, e.g., the 
CRNs given in mi and so it can make non-trivial predictions about the long¬ 
term behavior of those CRNs. In particular, the CRN N of Figure [T] satisfies 
the sufficient condition. Moreover, this result can also be used as a tool for 
engineering CRNs that perform deterministic computations (independent of the 
probabilities), such as in the computational model of [ 4 ]. Indeed, such CRNs 
generally require relatively simple long-term behavior which may be partially 
verified by Theorem [J. 

Theorem [J is inspired by the main technical result of [T] (which in turn was 
inspired by the main result of m). which provides another sufficient condition 
for the non-applicability of non-terminal reactions for recurrent configurations. 
However, there are a number of differences between both results. First, Theo- 
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rem [l| is derived in a basic combinatorial setting using notions from Petri net 
theory such as the notion of T-invariant, without considering stochastics. In 
contrast, the intricate proof of the main result of [ 1 ] is derived in a very differ¬ 
ent setting that uses non-trivial arguments from both mass action kinetics and 
stochastics. Secondly, we show examples where the main result of [1] is silent, 
while Theorem [l| makes a prediction. In fact, we conjecture that the main re¬ 
sult of □ is a special case of Theorem [l|. We compare both results in detail in 
Section y. While we focus in this paper on recurrent configurations of CRNs, we 
mention that the related concept of recurrent CRN has been investigated in [T3] . 

Deficiency is a useful and well-studied notion to classify CRNs. With Theo¬ 
rem [l| in place we consider at the end of Section 0 (and similar as done in [T]) 
its consequences for the case of CRNs with deficiency one (cf. Corollary!^. 

While formulated in terms of CRNs, the results in this paper equally apply to 
Petri nets, which is a very well studied model of parallel computation, see, e.g., 
M- Using the “dictionary” provided for the reader with a Petri net background 
(see Subsection l2.2l) . it is straightforward to reformulate the results in this paper 
in terms of Petri nets. 


2 Standard graph and CRN/Petri net notions 

2.1 Preliminaries 

Let N = {0,1,...}. Let X and Y be arbitrary sets. The set of vectors indexed 
by X with entries in Y (i.e., the set of functions : X T) is denoted by Y^. 
For v,w G N^, we write < w if v{x) < w{x) for all x € X. Moreover, we write 
V < w if V < w and v ^ w. The support of v, denoted by supp(z;), is the set 
{x G X \ v(x) > 0}. For finite sets X and Y, a X x Y matrix A is a matrix 
where the rows and columns are indexed by X and Y, respectively. 

We consider digraphs G = (V,E,F) where V and E are finite sets of vertices 
and edges and F : E ^ V'^ assigns to each edge e G E an ordered vertex pair 
(u,n). We denote V by V(G) and E by E{G). The incidence matrix of G is the 
V{G) X E{G) matrix A where for e G E with F{e) = {v,w) we have entries 
A(v,e) = —1, A{w,e) = 1, and A(u,e) = 0 for all m G U \ {w,w} if n ^ w, 
and A[u, e) = 0 for all u G V if v = w. The number of connected components 
of a digraph G is denoted by c(G). It is well known that the rank r{A) of the 
incidence matrix A of a digraph G is equal to \V\ — c{G) (where it does not 
matter over which field the rank is computed m Proposition 5.1.2]). From now 
on we let the field Q of rational numbers be the field in which we compute. 

A walk TT in G is described by (particular) strings over E. Let ^(tt) denote the 
Parikh image of tt, i.e., G where {<P{T)){e) is the number of occurrences 
of e in TT. We write supp(7r) = supp(^(7r)), i.e., supp(7r) is the set of elements 
that occur in tt. The vectors v of ker(A) flN^ describe the cycles of G, i.e., they 
describe the Parikh images of closed walks in G. 

For convenience we identify a digraph G with its V{G) x E{G) incidence 
matrix. Hence, we may for example speak of the rank r(G) of G. We say that 
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e G E{G) is a bridge if e is not contained in any closed walk of G. The induced 
subgraph G' of G with respect to X C V{G) is the digraph G' = {X,E',F') 
where E' is the preimage of under F and F' is the restriction of F to E'. 
A strongly connected component {SCC, for short) is an induced subgraph G' of 
G with respect to X C V(G) such that G" contains no bridge and X is largest 
(with respect to inclusion) with this property. 

2.2 CRNs and Petri nets 

We now recall the notion of a chemical reaction network. 

Definition 1. A chemical reaction network (or CRN for short) N is a 3-tuple 
(S, R, F) where S and R are finite sets and F is a function that assigns to each 
r £ R an ordered pair F{r) = {v,w) where v,w £ N'®. Vector v is denoted by 
in(r) and w by out(r). 

The elements of S are called the species of N, the elements of R are called the 
reactions of N, and F is called the reaction function. For a reaction r, in(r) and 
out(r) are called the reactant vector and product vector of r, respectively. 

It is common in the literature of CRNs to omit the function F and have R 
as a set of tuples {v,w). However, this would not allow two different reactions 
to have the same reactant and product vectors (such situations are common in 
Petri net theory). 

In CRN theory, it is common to write vectors in additive notation, so, e.g., 
if 5 = {A,B,G}, then A + 2B denotes the vector v with z;(A) = 1, v{B) = 2, 
and v{G) = 0. 

Example 1. Consider the CRN N = {S,R,F) with S = {A,B}, R = {a,b}, 
F{a) = (A + B,2B) and F{b) = {B,A). This CRN is taken from [TS] (see also 
m)- This example is the running example of this section. 

We now define a natural digraph for a CRN N, called the reaction graph of 
N. The name is from [10], and the concept is originally defined in |7]. 

Definition 2. Let N = {S,R,F) be a CRN. The reaction graph of N, denoted 
by TZn, is the labeled digraph (y,R,F) with V = {in(r) | r G i?} U {out(r) | r G 
R}. 

Note that in the reaction graph each reactant and product vector becomes a 
single vertex. The vertices of the reaction graph are called complexes. The re¬ 
action graph of the CRN N of our running example (Example [ij) is depicted in 
Figure yj. 


A + B 2B B —^ A 


Fig. 3. The reaction graph of the CRN of Example Q. 
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A configuration c of -/V is a vector c G N'®. Let r € R. We say that r can fire 
on c if in(r) < c. In this case we also write c — c' where c' = c — in(r) + out(r). 
Note that c' is a configuration as well. Moreover, we write c —)• c' if c — c' for 
some r £ R. For t £ R* (as usual, R* is Kleene star on R) we write c —c' if 
c Cl • • • — c' where r = n • • • and £ i? for alH £ n}. The 

reflexive and transitive closure of the relation —>■ is denoted by —If c c', 
then we say that c' is reachable from c. We say that a configuration c is recurrent 
if for all c' with c —>■* c' we have c' —>■* c. Note that if c is recurrent and c —>■* c', 
then c' is recurrent. 

Example 2. Consider again the running example. We have, e.g., 2A + B — 

2A + B. However, 2A + B is not recurrent as 2A + B ZA and in configuration 
ZA no reaction can fire. In fact, the recurrent configurations of N are precisely 
those that do not contain any B. Indeed, assume c is recurrent. Then we can fire 
b until we obtain a configuration c' that does not contain any B. No reaction 
can fire for c' and so c = c' since c is recurrent. 

The definition of a CRN is equivalent to that of a Petri net [14]. In a Petri net, 
species are called places p, reactions are called transitions^ and configurations 
are called markings. A Petri net is often depicted as a graph with two types of 
vertices, one type for the places and one for the transitions. The Petri net-style 
depiction of the running example is given in Figure 0. The round vertices are 
the places and the rectangular vertices are the transitions. We use in this paper 
several standard Petri net notions, which are recalled in the next subsection. 



Fig. 4. The Petri net-style depiction of the running example. 


2.3 P/T-invariants 

The notions of this subsection are all taken from Petri net theory [14]. We first 
recall the notion of an incidence matrix of a CRN, which is not to be confused 
with the notion of an incidence matrix of a digraph (as recalled above). In fact, 
we will compare in the next subsection the incidence matrix of a CRN with the 
incidence matrix of its reaction graph. 
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Definition 3. For a CRN N = {S,R,F), the incidence matrix of N, denoted 
by In , is the S x R matrix A where for each r G R the column of A belonging 
to r is equal to out(r) — in(r). 

Example 3. Consider again the CRN N of the running example. Then 


a b 



Note that if c c', then c' = c + 2jv^(t), where 'P{t) denotes again the 
Parikh image of r. 

A u e N® is called a P-invariant of N if v'^In = 0 (here 0 denotes a zero 
vector of suitable dimension indexed by R). Similarly, v G is called a T- 
invariant of N if Inv = 0, i.e., v G ker(Xjv)0 A P-invariant or T-invariant are 
also sometimes called P-semiflow and T-semiflow, respectively, in the literature. 
Observe that if c c', then I’{t) is a T-invariant if and only if c' = c. A CRN 
N is called conservative if there is a P-invariant v such that supp(u) = S. Also, 
N is called consistent if there is a T-invariant v such that supp(u) = R. 

A CRN N is said to be structurally bounded when for every configuration c, 
there is a fcc G N such that for each configuration c' with c —>■* c' we have that 
each entry of c' is at most kc- Note that for a structurally-bounded CRN, the 
number of different configurations reachable from a given configuration is finite, 
and so for each configuration c, there is a recurrent configuration reachable from 
c. In this way, one often informally views the recurrent configurations as the 
possible states of the CRN in “the long term”. 

The following result is well known, but for completeness we recall its short 
proof. 

Proposition 1 (|Tl]). Let N be a CRN. If N is conservative, then N is struc¬ 
turally bounded. 

Proof. Let v G be a P-invariant with supp(u) = S and let c be a configuration. 
Let c c' for some r G R*. We have c' = c -I- In‘L{t). Thus v'^c' = v'^c -\- 
v'^lNd>{T) = v'^c and so for all s G S, v{s)c'{s) < v'^c and therefore c'(s) < 
v’^c/v{s). □ 

Example 4- The CRN N of the running example is both conservative and consis¬ 
tent. Indeed, any v G with v(A) = v{B) > 1 is a P-invariant with supp(u) = S 
and any w G with v{a) = v{b) > 1 is T-invariant with supp(u) = R. 

2.4 Deficiency 

The notions that we recall in this subsection are originally from chemical reaction 
theory (and are less studied within Petri net theory). 

^ The P and T in P/T-invariant are short for Place and Transition (from Petri net the¬ 
ory) . We choose to use these well-known names instead of calling them “S-invariant” 
and “R-invariant” for Species and Reaction, respectively. 
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Let N = {S, R, F) be a CRN and let V = {in(r) | r G i?} U {out(r) | r G i?}. 
We denote by the S xV matrix with for all s G S' and v gV, entry >’Ar(s, v) 
is equal to v{s). 

The next lemma relates the incidence matrix In of a CRN N with the 
incidence matrix of the reaction graph TZn of N. 

Lemma 1 (Section 6 of [8]). Let N = {S,R,F) be a CRN. Then In = 
yNR-N- 

In the above equality, TZn denotes the incidence matrix TZn and not the graph. 

Proof. Let V = {in(r) | r G i?} U {out(r) | r G R}. Let p G P and r G R. 
Then lN{p,r) = (out(r) - in(r))(p) = 3^w(p, out(r)) • 1 + 3^7v(p, in(r)) • (-1) = 
Y.x(ivyN{P^x)nN{x,r) =yN'R-N- □ 

As a corollary to Lemma [H, we have the following. 

Corollary 1 I[10jl. Let N = {S,R,F) be a CRN. Then ker(7?.Ar) C ker(lAr). 

The vectors v of ker(7?.Ar) Cl N^, which are T-invariants by Corollary [l|, are 
called closed T-invariants [3]- Recall that the vectors v of ker(72.Ar) flN^ describe 
the cycles of TZni and so for each closed T-invariant v of N, supp(?;) does not 
contain any bridge of TZn- Since each of the entries of a T-invariant is nonnega¬ 
tive, the linear space ker(Iw) does not necessarily have a basis consisting of only 
T-invariants, see Example @ below. 

The deficiency S{N) of a CRN N is r(7?.jv) — r(Ijv)- By Corollary [l|, S{N) 
is non-negative. Thus, one may view S{N) as a measure of the difference in 
dimensions between ker( 7 ^ 7 v) and ker(l 7 v)- The former is determined only by 
the structure of the reaction graph (ignoring the identity of the vertices), while 
the latter also incorporates the relations that rely on the identities of the vertices 
of the reaction graph. 

Recall from Subsection 12 .II that r{TZN) = \V{TZn)\ —c{TZn)- Hence, we have 
6{N) = \V(TZn)\ — c{TZn) — r(lN) |9I7) . Note that if d{N) = 0, then every 
T-invariant of N is closed and ker(7^Ar) = ker(l 7 v)- 


A + B 2B A —^ B 


Fig. 5. The reaction graph of a CRN discussed in Example 


Example 5. In the running example, ker(77.jv) only contains the zero vector, 
while ker(lAr) contains all scalar multiples of the vector w with w{a) = w{b) = 1 . 
Thus ker(lAr) has a basis consisting of only T-invariants. Moreover, &{N) = 1. 
Alternatively, the reaction graph TZn has 4 vertices and 2 connected components 
and r(lAr) = 1. Thus, 5{N) = 4 - 2 - 1 = I. 

If we consider the CRN N' of Figure d, then ker(72.Ar/) also only contains 
the zero vector, while ker(Iw') contains all scalar multiples of the vector w with 
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w{a) = —w{b) = 1. Again, S{N') = 1, however the only T-invariant of ker(Ijv/) 
is the zero vector. 

3 Dominance and non-closed T-invariants 

Note that there is a natural partial order for the set of SCCs of a graph: for 
SCCs A and Y, we have A ^ F if there is a path from a vertex of F to a vertex 
of A. We now consider a different partial order, denoted by <d, for the SCCs of 
a reaction graph of a CRN. 

Let A be a CRN. For SCCs A and F of TZn we write A <d F if there are 
vertices a: of A and y of F such that x <y. 

Lemma 2. Let N = {S,R,F) be a structurally-bounded CRN. Then the <d 
relation between SCCs of TZn is a partial order. 

Proof. The <d relation is obviously reflexive and transitive. To show that <d is 
antisymmetric, let X <dY and Y <d X for some SCCs A and F oUZn- Hence 
there are vertices xi and X 2 of A and yi and j /2 of F such that xi < yi and 
t /2 F X 2 . Let TTi be a path from xi to X 2 and let ^2 be a path from 1/2 to yi 
in TZn- Then xi + 7/2 X 2 + j /2 X 2 + yi- If A F, then xi < yi and 
y 2 < X 2 . Thus we have xi + 7/2 < X 2 -\-yi, and so N is not structurally bounded 
— a contradiction. □ 

For SCCs A and F we write A <d F if A F and A 7 ^ F. We say that A 
dominates Y when A <d Y. For a set S of SCCs, we let min<^( 45 ) C 5 be the 
set of elements of S that are minimal with respect to the <d relation among all 
the elements of S. 

Let us define for a SCC A of TZn, out(A) = {r G E{TZn) \ in(r) G 
F(A),out(r) ^ F(A)}. We call A terminal if out(A) = 0. We call a reaction 
r (complex x, resp.) terminal if r G E{X) {x G F(A), resp.) for some terminal 
SCC A of 7^Ar. 

We will consider the minimal set X of non-terminal SCCs that dominates all 
other non-terminal SCCs. In other words, if we let JV be the set of non-terminal 
SCCs, then X = min<^(A/"). 

Let B be the set of bridges oi TZn- The exit set of a set S of non-terminal 
SCCs, is a set F C i? with both \Z\ = l^l and \Z n out(A)| = 1 for all A G 5. 
In other words, Z contains exactly one bridge of out(A) for each X € S- 

Assuming the existence of a non-terminal reaction that can fire for some 
recurrent configuration c, the main result of this paper ensures the existence of 
certain sequences t with c' -j-'^ d for some configuration c' reachable from c. For 
each exit set F, there exists such a t that avoids all bridges outside Z and, at 
the same time, uses the bridges of Z whenever possible. As a consequence, each 
of the sequences r corresponds to a T-invariant v = ^{t) that have zero entries 
for the bridges outside Z and nonzero entries for some of the bridges inside Z. 
We will show that for various CRNs this necessary condition allows one to show 
that only terminal reactions can fire for all its recurrent configurations. 
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The proof idea is the following. Let us start with a recurrent configuration c. 
While traversing the configuration space by applying reactions starting from c, 
we need never choose a bridge of 7 ?.at going out of a SCC X that is dominated 
by some Y (i.e., Y <d X). Indeed, if cc G V{X) and y G V{Y) with y < x, then 
we may walk inside X to x and y < x implies that any reaction r with in(r) = y 
can fire for x. In this way we also avoid taking a reaction r' with in(r') = x. 
Moreover, walking out of Y can be done by taking any of the bridges. We choose 
the one from the exit set Z. Now, eventually, our path inside IZn will lead to 
a terminal vertex. However, since c is recurrent, we can go back to c. If a non¬ 
terminal reaction can fire for c, then this means that we can iterate this process 
(walking along bridges, etc.). Structural boundedness finally ensures that the 
configuration space is finite and so, we must eventually repeat a configuration 
that closes the “circuit”. 

We are now ready to formulate the main result of this paper. 

Theorem 1. Let N = {S,R,F) be a structurally-bounded CRN, and let X = 
min<^(j\/’), where N is the set of non-terminal SCCs ofTZjsf. Let B be the set of 
bridges oflZN- Let L be the set of all non-terminal reactions r oflZ^ such that 
there is a non-terminal reaction r' ofTZ]\[ with in(r') < in(r). 

If some non-terminal reaction can fire for some recurrent configuration c, 
then for each exit set Z of X, there is a t G R* such that 

1. T contains no reactions from {B\Z)\J L, 

2 . T = TTiCTi • • • 7r„tT„ where each tt^ is a path in TZn from a non-terminal vertex 

to a terminal vertex and each ai is a sequence of terminal reactions, and 

3. c' c' for some recurrent configuration c' reachable from c. 

Proof. Assume that some non-terminal reaction ri can fire for some recurrent 
configuration c. Let Z C B he an. exit set of X. 

Let Y be the SCC containing vertex in(ri) and let si G i?* be a shortest 
path in Y from in(ri) to some vertex y with x < y for some vertex x of some 
X G X. Observe that supp(si) flL = 0 (if F G or ri G L, then si is the empty 
string). Let S 2 bx G i?* be a path in TZx such that S 2 is a path in X from x and 
bx G Z f] out(A). Then ai = riSiS 2 bx can fire for c. We repeat this procedure 
for vertex out(6x)- In this way, we obtain a path 0102 - By iteration, we obtain 
a path TTi = ai ■ ■ ■ ai from in(ri) to a terminal vertex x. Note that tti can fire 
for c. Since c is recurrent, there is a sequence cri of terminal reactions such that 
Ti = TTiCTi has the property that c — c' where some non-terminal reaction r 2 
can fire for c'. Note that both ri contains no reactions from (B \ Z) U L. Also, 
TTi is a path inIZx from a non-terminal vertex to a terminal vertex. 

We repeat the above described procedure (that constructed n) for con¬ 
figuration c' and vertex in(r 2 ), to obtain (by iteration) an infinite sequence 
T = T1T2 • ■ • . Since r is infinite and N is structurally bounded, there is a con¬ 
figuration Cr such that c —Cr ^"^ippp Cr and rioop = n - ■ - tj for some i < j. 
Note that, by the construction of r, rioop contains no reactions from {B\Z)UL, 
and so we are done. □ 

We illustrate Theorem [H through a couple of examples. 
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A — C—E + D—^A + C 


Fig. 6. The reaction graph of the CRN of Example 


Example 6. Consider the CRN N of FigureIt is easy to verify that c = A-\-C 
is a recurrent configuration. Moreover, there is a non-terminal reaction r that 
can hre for this configuration (take r = a or r = b). Note that there is only one 
exit set Z for X, which is Z = B = {a,b,d}. By Theorem [U there is a r £ i?* 
such that (1) r contains no reactions from {B\Z)UL, (2) r is a sequence of paths, 
each going to a terminal vertex, and (3) c' c' for some recurrent configuration 
c' reachable from c. Indeed, we can choose, e.g., r = abd and E = Ax C. 


A + C —^E + C E + D — ^-^A + D 


Fig. 7. The reaction graph of the CRN of Example 0. 


We now give another example. 

Example 7. Consider the CRN N of Figure 0. It is easy to verify that c = 
A+C+D is a recurrent configuration. Moreover, there is a non-terminal reaction 
r that can fire for this configuration (take r = a). We have that Z = B = {a, b} 
is the unique exit set Z ioi X. We notice that t = ab and c' = A + C + D satisfy 
the conditions of Theorem [TJ Indeed, we have t = Ticrir 2 cr 2 with ri = a, cti = e 
(the empty string), T 2 = &, and a 2 = e. Note that if N contained the additional 
reaction A + D E + D, then t = ab and c' = A + C + D would again satisfy 
the conditions of Theorem [T1 where r = riui with ti = a and ci = b. 

Considering the non-closed T-invariant v = d>{T) with r from Theorem^, we 
have the following corollary to Theorem [l|. Note that Condition [5] of Theorem 0 
implies that supp(t!) contains a bridge, and therefore v{z) fy 0 for some z € Z. 

Corollary 2. Let N, X, B, and L be as in Theorem\^. 

Assume there is an exit set Z of X such that there is no non-closed T- 
invariant v with (1) v{x) = 0 for all x € {B \ Z) U L and (2) v{z) fy 0 for some 
z G Z. 

Then no non-terminal reaction can fire for any recurrent configuration of N. 

We remark that, in view of Theorem 0 , Corollary 0 can be strengthened 
by replacing the condition v{z) fy 0 for some z G Z with the stronger (but 
more involved) condition that says that the (occurrences of the) non-terminal 
reactions of v form a set of paths where each path ends in a terminal vertex. 

Note that since closed T-invariants v cannot contain bridges, we may without 
loss of generality remove the condition that v is “non-closed” in Corollary 0 . 
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We use Corollary 0 to determine whether no non-terminal reaction can hre 
for any recurrent configuration of a CRN. While non-closed T-invariants have a 
central role in Corollary 0, curiously, this notion from [3] has been given only 
modest attention in both the Petri net theory and the CRN theory. 

For a given exit set Z oi X, one can verify using linear programming in 
polynomial time whether or not there is a non-closed T-invariant v with the 
properties of Corollary!^. While in general there may be an exponential number 
of exit sets (exponential in the number of reactions) to check, in many cases 
the number of exit sets is severely constraint and in these cases the sufficient 
condition of Corollary [3 is computationally efficient. 


A 




D D + E > J + H 


A + H A + E C + H —^C + E 


Fig. 8. The reaction graph of the CRN of Example 


We now give some examples to illustrate Corollary 

Example 8. Consider the CRN N of Figure 0. This CRN is a simplification of a 
CRN from biology studied in [15] (see also my We have 


a 

A /-I 


Ijv = 


J 

C 

D 

E 

H 


1 
0 
0 
0 
V 0 


-1 -1 


d 

0 

1 


-1 -1 


/ 

0 

1 

0 

-1 

-1 

1 


9 

0 

0 

0 

0 

1 

-1 


h 

0 

0 

0 

0 

1 

-1 


\ 


/ 


It is easy to verify that the sum of the rows of Ijv is the zero vector and so N is 
conservative. Consequently, N is structurally bounded. It turns out that ker(l 7 v) 
is of dimension 4 and is spanned by T-invariants. In fact, one can verify that 
ker(Ijv) is spanned by the two closed T-invariants wi = <l>{ab) and W 2 = d){cd) 
together with the two non-closed T-invariants vi = <P{gfce) and V 2 = d>(hfce). 
We remark that A + H + D -^afce + D axvdCAH + D C + H + D. 

Thus S{N) = 2. Note that B — {e,f,g,h} is the set of bridges of TZn- Let X 
be the set of non-terminal SCCs of TZn that are minimal with respect to <d- 
We notice that Z = {e, /} is the only exit set of X. Also L = { 5 , h}. Now, 
the non-closed T-invariants vi and V 2 are witnesses that there is no non-closed 
T-invariant v with both (1) v{g) = v(h) = 0 (Aote that {B\Z)\JL = {g, h}) and 
(2) either v(e) or v{f) nonzero. By Corollaryl2|, for every recurrent configuration 
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no non-terminal reaction can fire. Since every reaction is non-terminal, for every 
recurrent configuration no reaction can fire. 


A + B -2-^ 2B B + C —^ A + C 


Fig. 9. The reaction graph of the CRN of Example [^. 


The next example shows that the converse of Theorem [l| does not hold. 

Example 9. Consider the CRN N of Figure]^. We show that no reaction can fire 
for any recurrent configuration of N. Let c be a recurrent configuration. If c does 
not contain any (7, then we can fire reaction a until we obtain a configuration 
d for which no more reactions can fire. Since c is recurrent, c = c' and we are 
done. If c contains at least one C, then we can apply reaction b until we obtain a 
configuration c" with only A’s and C’s. Hence no reaction can fire for c". Since 
c is recurrent, we have c = c" and we are done. 

However, for c = A -|- H -|- C we have c —c with r = ab. We notice that 
Z = {a, b} is the only exit set of X and {B\Z)\JL = 0. Thus r trivially contains 
no reactions from {B\Z)\J L and t = tti'K2 with paths tti = a and 7r2 = 6 in 
from non-terminal vertices to terminal vertices. This shows that the converse of 
Theorem [l| does not hold. 

We remark that if we remove species C from reaction 6, then Corollary 
(and Theorem [J) would have been applicable to show that no (non-terminal) 
reaction can fire for any recurrent configuration of N. 

We now consider the case where the deficiency is 1. This severely restricts 
the structure of the non-closed T-invariants. 

Lemma 3. Let N = (S', i?, F) be a consistent CRN with 6{N) = 1. Then for all 
non-closed T-invariants v, supp(u) contains every bridge ofTZi\j- 

Proof. Let u be a non-closed T-invariant and let 5 be a bridge of 7 ^at. Since 
N is consistent, there is a T-invariant w with supp(w) = R. Since b G R, w is 
non-closed. Thus v,w £ ker(XAr) \ker(7?.Ar). Since S{N) = 1 and supp(z) cannot 
contain any bridge for z £ ker(7?.Ar), we have that b £ supp(w) if and only if 
b £ supp(u). Hence b £ supp(u). □ 

The next result is essentially Theorem 3.5 of the supplementary material of 
[T] (although there it is stated in terms of notions from mass-action kinetics and 
stochastics), and follows directly from Theorem [l| and Lemma 0. 

Corollary 3 ([T]I. Let N be a structurally-bounded and consistent CRN with 
S{N) = 1. If there are non-terminal vertices x and y such that x < y, then for 
all recurrent configurations c, none of the non-terminal reactions can fire. 
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Proof. Assume there are non-terminal reactions x and y such that x < y and 
assume to the contrary that some non-terminal reaction r can fire for some 
recurrent configuration c. By Lemma0, for all non-closed T-invariants u, supp(u) 
contains every bridge of IZn- Hence, by Theorem [l|, every non-terminal SCC of 
TZn is minimal with respect to <d among the non-terminal SCCs of TZn — a 
contradiction by the existence of x and y (note that x and y cannot be vertices 
of the same SCC since N is structurally bounded). □ 

Example 10. Consider the CRN N of the running example of Section [3. Recall 
that N is conservative, and therefore N is structurally bounded. Also recall that 
5{N) = 1. By Corollary 0, no non-terminal reaction can fire for any recurrent 
configuration c of N. Since all reactions of N are non-terminal, no reaction can 
fire for any recurrent configuration c of A^. 


4 Using rates 

This paper is inspired by the main technical result of [T] (cf. Theorem 3.3 of the 
supplementary material of m)- In this section we recall its result. First we recall 
a particular matrix. Let R>o (R> 0 ; resp.) be the set of nonnegative (positive, 
resp.) real numbers. 

Definition 4. Let N = {S, R, F) be a CRN. Let V = VifR-j^i) and let k € R> 0 ' 
We denote by ICn.k, the S' x F matrix where for each x £V the column of 
belonging to x is equal to in(r)=x ' (out(r) — in(r)). 

The value K{r) in Theorem@may be interpreted as the “rate” of reaction r. Note 
that the definition of is closely related to the definition of (Definition^). 
We are now ready to formulate the main technical result of [T]. 

Theorem 2 ([!]). Let N = (S,i?,F) be a conservative CRN and V = V(TZn). 
Let L be the set of non-terminal vertices v of TZn such that there is a non¬ 
terminal vertex v' oflZ^ with v' < v. Assume that L 7 ^ 0. 

If some non-terminal reaction can fire for some recurrent configuration c, 
then for all k there is a w G ker(/CAr_„) fl R^o with supp(w) n L = 0 and 

there is a non-terminal vertex x with x G supp(w). 

Theorem 0 is proved in [T] using both intricate probabilistic arguments and 
methods from mass action kinetics. In [T], the theorem is unnecessarily stated 
in a probabilistic fashion using the notion of “positive recurrent configuration” 
for stochastically modeled CRNs: it can be stated in a deterministic way (see 
Theorem above) by realizing that the configuration space is finite for a given 
initial configuration in a structurally-bounded CRN. This deterministic formu¬ 
lation and the discrete model (in contrast to mass action) triggered the search 
of this paper for a combinatorial explanation of this result. We invite the reader 
to compare the proof techniques used to prove Theorem in [T] and Theorem [l| 
in this paper. In [1], Corollary H is proved using Theorem @ while in this paper 
it is shown using Theorem [l|. 
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2B A + B 2A 


Fig. 10. The reaction graph of the CRN of Example [ill. 


Note that if L = 0, then Theorem @ is silent. We now show an example with 
L = 0 where Corollary can be applied. 

Example 11. Consider the CRN N of Figurefl^ Note that N is conservative with 
w{A) = w{B) = 1 as a witness. The only T-invariants v oi N are those where 
v{a) = v{b). Let Z = {a} be an exit set of X. Then there is no non-closed T- 
invariant v with v{h) = 0 and v(a) ^ 0. By Corollary^, no non-terminal reaction 
can fire for any recurrent configuration c of Since all reactions of N are non¬ 
terminal, no reaction can fire for any recurrent configuration c of N. Indeed, 
one observes that the recurrent configurations of N are those configurations 
containing either only A’s or only _B’s, for which a and b cannot fire. 

We conjecture that the assumption L 7 ^ 0 can be removed from Theorem 0. 
In case L 7 ^ 0 is removed from Theorem 0, then Theorem 0 also predicts that 
no non-terminal reaction can fire for any recurrent configuration of the CRN of 
Example El Next, we give an example with L 7 ^ 0, where Corollary @ can be 
applied but Theoremis silent. 


a 

A + D' ^ B + D—+ D 

~ K — 

b 


B + E' A + E—-^C + E 

~ K — 

e 


2A + D — -^ 3F 


Fig. 11. The reaction graph of the CRN of Example [l^ . 


Example 12. Consider the CRN N of Figure I 111 Note that N is conservative 

+ D < 2A + D and 


with 

w{X) 

= 1 for all 

species X as 

a witness. 

N^ote that 

so L 

7 ^ 0 in 

Theorem y. Let k € R>o- We have 

Xn,k = 



A + D 

B + D 

B + E 

A + E 

2A + D 

A 

/ 

—K(a) 

K{b) 

K{d) 

-K(e) - k(/) 

- 2 K(g) 

B 


K(a) 

— k(&) — K 

(c) -K{d) 

«(e) 

0 

C 


0 

k{c) 

0 

k(/) 

0 

D 
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0 

0 

0 

-Ea) 

E 


0 

0 

0 

0 

0 

E 

1 

0 

0 

0 

0 

3K{g) 


0 \ 
0 
0 
0 
0 
0 y 


Let w S ]R>q with K{a)w{A-\- D) = K{d)w{B + E) > 0 and w[x) = 0 for all other 
X & V . Then w S ker{ICN,it) C with X € supp(?u) for some non-terminal 
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vertex x and supp(r(;) H L = 0. Thus Theorem [3 is silent. On the other hand, 
none of the non-closed T-invariants of N contains a bridge and so by Corollary^, 
no non-terminal reaction can fire for any recurrent configuration of N. 

Conversely, despite trying numerous examples, we could not find an example 
where Theorem^ predicts that no non-terminal reaction can fire for any recurrent 
configuration, but where Theorem [l| is silent. 

5 Discussion 

Based on structural properties of CRNs, the main result of this paper (cf. The¬ 
orem [l|) provides a sufficient condition to analyze the long-term behavior of 
CRNs. While its proof is using basic combinatorial arguments, the result is pow¬ 
erful enough to apply to a large class of CRNs. Also, the sufficient condition is 
computationally-efhcient to verify for many CRNs. Another such sufficient con¬ 
dition is shown in [1], cf. Theorem [H. We have shown examples of CRNs where 
Theorem [H is applicable while Theorem 0 is silent. 

Given that discrete CRNs are equivalent to Petri nets, it is curious that the 
corresponding research areas of CRN theory and Petri net theory have evolved 
almost independently. In this paper we shown that notions from Petri net theory 
(in particular, T-invariance) are useful for CRN theory. Similarly, notion such 
as dehciency, originating from CRN theory, are useful for Petri net theory. At 
the interface of these two notions is the scarcely-studied notion of non-closed 
T-invariant, which is crucial in the sufficient condition of Corollary 0. This illus¬ 
trates that both research areas can significantly profit from each other. 

An open problem is resolving whether Theorem is indeed a special case 
of Theorem [l|. Another open problem is to somehow strengthen Theorem [l| in 
a natural way to make it applicable for CRNs such as the one presented in 
Example 0. 

A further research direction is to incorporate probabilities. One may asso¬ 
ciate a probability to each T-invariant by multiplying the probabilities of the 
corresponding reactions. An open problem is to find a probabilistic version of 
Theorem [l| to make predictions about long-term behavior of probabilistic com¬ 
putational models of CRNs, such as the models of mm- 
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